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Considering systems of self-propelled polar particles with nematic interactions (“rods”), we com¬ 
pare the continuum equations describing the evolution of polar and nematic order parameters, 
derived either from Smoluchowski or Boltzmann equations. Our main goal is to understand the dis¬ 
crepancies between the continuum equations obtained so far in both frameworks. We first show that 
in the simple case of point-like particles with only alignment interactions, the continuum equations 
obtained have the same structure in both cases. We further study, in the Smoluchowski framework, 
the case where an interaction force is added on top of the aligning torque. This clarifies the origin 
of the additional terms obtained in previous works. Our observations lead us to emphasize the need 
for a more involved closure scheme than the standard normal form of the distribution when dealing 
with active systems. 


I. INTRODUCTION 

“Self-propelled rods”, i.e. elongated objects spending 
energy to displace themselves typically along their long 
axis, are among the most generic and ubiquitous objects 
studied in active matter physics. Living and inert exam¬ 
ples abound: elongated bacteria swimming or crawling 
on a surface [l|, chemically-propelled micro- and nano¬ 
rods 0, biofilaments displaced by molecular motors [1- 
0, shaken granular particles Q, etc. Their main inter¬ 
action often amounts to alignment due, e.g., to direct 
contact via collisions. In many cases the fluid in which 
the self-propelled rods move can be treated as inert, sim¬ 
ply providing friction. This “dry active matter” settiM 
nevertheless shows non-trivial collective properties 
[Tsjl . Among these, the emergence of dense (almost close 
packed), polarly-oriented clusters has been noted repeat¬ 
edly and is a quite natural outcome of aligning collisions: 
when two rods meet, they often either end up aligned 
or anti-aligned (“nematic” alignment). In the first case, 
they will stay alongside each other for quite a while, even 
in the absence of attractive interactions, and may recruit 
others, forming clusters. 

Numerical work on various models of moving elon¬ 
gated objects interacting via steric repulsion has revealed 
the possibility of complex (and not yet fully understood) 
phase diagrams, even for simple rigid rods. The ubiquity 
of dense clusters, forming even at low global densities, 
seems a priori to be a major difficulty in building theoret¬ 
ical, continuum descriptions of these systems since most 
approaches so far rely on the hypotheses that only bi¬ 
nary interactions need to be considered and/or that rods 
decorrelate between collision events. Nevertheless, such 
continuum descriptions have been considered, notably by 
us, in the past. 

Theoretical approaches typically yield continuum (or 
“hydrodynamic”) equations governing a nematic and a 


polar order field, as well as a continuity equation for the 
density field. By ’continuum equations’, we mean here 
a reduced set of evolution equations for conserved and 
order parameter fields. The first derivations of such sets 
of equations were performed by Baskaran and Marchetti, 
who treated explicit collisions between thin rods dill, 
while a different route, treating alignir^ point-like par¬ 
ticles, was followed by Peshkov et al. [^. These works 
predicted the emergence of global nematic order at the 
deterministic level, that is, without including noise in the 
continuum equations. 

Baskaran and Marchetti first studied the case of in¬ 
teracting self-propelled rods assumin g ov erdamped mi¬ 
croscopic dynamics from the outset d- ^ second 
paper they included inertia in the microscopic dynamics 
[l7|, with linear and angular momentum transfer during 
collisions between rods. They derived a Fokker-Planck 
equation for the joint probability distribution of posi¬ 
tions and velocity, and then took the overdamped limit 
of the kinetic equation, obtaining a Smoluchowski equa¬ 
tion for the one-particle probability density. Continuum 
equations for the density, polar and nematic fields were 
then obtained from the Smoluchowski equation, and the 
linear stability of the basic homogeneous states was stud¬ 
ied. This derivation yields important modifications of 
the Smoluchowski equation that result in nonlinearities 
in the continuum equations that are not obtained when 
considering overdamped dynamics from the outset. Here 
we focus on the equations reported in m. 

Following the kinetic approach pioneered by [l^ for 
polar constant-speed point-like particles aligning ferro- 
magnetically, Peshkov et al. treated the case where the 
same particles align nematically, i.e. anti-align when their 
incoming angle is larger than ^ [2^ . They obtained well- 
behaved nonlinear partial differential equations whose so¬ 
lutions were shown to be in good qualitative agreement 
with the Vicsek-like model they were derived from. 
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The sets of equations resulting from these works bear 
strong similarities, but also differences. In this paper, we 
explain and discuss the origin of these differences. 


II. MODELS AND MAIN GOALS 

We consider self-propelled polar particles moving in 
two dimensions at constant speed vq. When isolated, 
their polarity angle 6 diffuses at long times. They align 
nematically with neighbors. To match previous studies 
[la flTl [^ . , we consider either pure diffusion of the 

angle 9 with rotational diffusion coefficient Dji, or ’run- 
and-tumble’ type dynamics where 9 is changed, with a 
probability A/j per unit time, into 9' = 9 + rj, where 
ry is a random variable of distribution P{ri). Moreover, 
we neglect positional diffusion in sections HIl and Hill and 
later reintroduce it in section HVl 

In the following, we compare continuum equations ob¬ 
tained starting from a Smoluchowski equation and from 
a Boltzmann equation. Both these kinetic equations gov¬ 
ern the time evolution of the probability /(r, 9, t) to have 
a particle at position r with polarity 9, at time t. Al¬ 
though derived from different starting points, these equa¬ 
tions bear strong formal similarities, and one may thus 
expect the resulting continuum equations describing the 
evolution of the relevant order parameter fields (here, the 
polar and nematic fields) to be similar, with changes af¬ 
fecting only the precise values of the coefficients. It was 
found, however, that the continuum equations obtained 
in |l6l . [l7| startin g fr om the Smoluchowski equation on 
one side, and in |20l . starting from the Boltzmann 
equation on the other side, possess many differences, as 
described below. The goal of this paper is on the one 
hand to emphasize the formal similarities of the two ap¬ 
proaches, by formulating them in a common framework, 
and on the other hand to outline the origins of the dif¬ 
ferences in the resulting continuum equations. 


A. Comparing Smoluchowski and Boltzmann 
equations 


side both equations contain a convective mass flux due 
to self-propulsion at speed uq, with e{9) = (cos0, sin0) 
a unit vector along the direction of self-propulsion. The 
speed Vq was set equal to 1 in [2^ (without loss of gen¬ 
erality), but is retained here for clarity. 

The content of the right hand side of Eq. 0 can be 
described as follows: Dji is a rotational diffusion coeffi¬ 
cient representing a Gaussian white noise process that re¬ 
orients the particle’s self-propulsion velocity at each time 
step; r is the mean-field torque exerted by the other rods, 
given by (to first order in spatial gradients) 

r(r,0,f) = J d9'Ki {9 - 9') f {r,9' ,t) 

+ Jd9'K2{9',9)-Vf{r,9',t), (3) 

where Ki and K 2 are collision kernels describing hard 
rods collisions. The second term in Eq. @ takes into 
account the difference in position of the center of the 
colliding rods due to their finite size (see Appendix . 
Finally, the mean-field force F is given by (to lowest order 
in gradients) 


F= [\9'G{9,9')f{r,9',t) , (4) 

J —TT 


where G is again a hard rod collision kernel describing 
linear momentum transfer in a collision. 

The content of the right hand side of Eq. m can be 
described as follows: /dif[/] is a generator of rotational 
reorientation through an arbitrary stochastic process ry. 


/dif[/] = A« [ d9'f{9')[{S2Ad'+ V-0))^ - 62 .( 0 '-0)] 

J —TT 

( 5 ) 

(with Afl a frequency of “tumbling” events) and /coi [/] is 
the collision integral defined as 


IcoAf]= ! d9i [ d92KB(92-9i)f(r,9i)f(r,92) 

J —TT J —TT 

X {62.{^(9i, 92) + rj — 9'))n — 52.(9 i — 9) (6) 


The Smoluchowski and Boltzmann equations for the 
probability density f(r,9,t) are respectively given by 

dtf + voe(9) ■ V/ = Dad^f - 9,(/r) - V ■ (fF) (1) 

dtf + voe(9) ■ V/ = /dif[/] + /coi[/] (2) 

(see below for notations). The Smoluchowski equa¬ 
tion 0 was derived in |l8l | by coarse-graining a micro¬ 
scopic model of long, hard thin rods undergoing colli¬ 
sions. That work also included translational diffusion, 
which is neglected here to simplify the comparison be¬ 
tween the two approaches (the effect of translational dif¬ 
fusion is discussed in Sec. lYl). The Boltzmann equa¬ 
tion @ was obtained in [ 2 ^ for point particles with pre¬ 
scribed nematic alignment interactions. On the left hand 


where 62 . is a generalized Dirac distribution taking into 
account the 27r-periodicity of angles. The notation (... )r) 
indicates an average over the noise distribution P(rj), and 
'I'(^ii ^ 2 ) is the direction of motion of particle 1 after the 
collision between particles 1 and 2, up to an additive 
noise ry (see Sect. IIIIl for details). We assume here that 
'I'(^ 11 ^ 2 ) favors nematic alignment. 

There are three differences between these two kinetic 
equations, (i) First, the reorientation events in the Boltz¬ 
mann approach are described by an arbitrary stochastic 
process with probability distribution P (rj). For a uni¬ 
form distribution P(r]) over [—7r,7r], one recovers run- 
and-tumble dynamics, while a (wrapped) Gaussian P (rj) 
yields Brownian rotational diffusion at rate A/j. In the 
Smoluchowski equation, in contrast, one has assumed 
Gaussian rotational noise from the outset. This slight 
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difference, however, has no influence on the structure of 
the resulting continuum equations. Moreover, when a 
van Kampen expansion is carried out in Eq. in the 
case of a distribution P{t]) with a small variance tr^ ^ 1, 
one finds /dif[/] — Dudgf, with Du = Ai^cr^, so that 
one recovers a rotational diffusion term as in Eq. O- 
(ii) Secondly, the Smoluchowski equation contains spatial 
gradients in both force and torque, while the collision in¬ 
tegral is local in the Boltzmann equation. This difference 
arises because the Smoluchowski equation describes rods 
of finite length I, necessitating an additional gradient ex¬ 
pansion in the size of the particles, while the Boltzmann 
equation assumes point particles, (iii) Finally, a third 
subtle distinction exists between the two theories and 
lies in the details of the collision kernels whose explicit 
forms are discussed below. In the Boltzmann description, 
combining the prescribed nematic alignment rule with 
the kinetics of collision leads to a kernel of mixed polar 
and nematic symmetry. On the other hand the collision 
kernel Ki associated to the torque in the Smoluchowski 
equation, calculated from hard rod collisions , has ne¬ 
matic symmetry in the limit of infinitely thin rods. As 
we shall show below, the nematic symmetry is also re¬ 
covered in the Boltzmann case by considering the limit 
of infinitely thin rods (while the opposite limit of quasi¬ 
circular particles was originally considered in HI). 


B. Comparing the derived macroscopic 
hydrodynamics 

Macroscopic continuum equations have been derived 
from both the above Smoluchowski and Boltzmann equa¬ 
tions in [la and in Hsini. A first step consists in trans¬ 
forming these equations into a hierarchy of field equations 
by introducing the Fourier expansion of /(r, 0, t) in 9: 

fk{r,t)= f (7a) 

J —TT 

OO 

fir,9,t) = — (7b) 

k— — oo 


comparison|3 

9tP + z;oRe(V*/i) = 0, (9a) 

dtfi = -Dnfi + Cs/r/2 - y (Vp + VV 2 ) 

- A'(/rV/i + /iVVi - /iV/r) , (9b) 

dtf2 = Ps/2 - y V/i - ^/rV/2 - ^/2V/r 

- ^/ 2 V*/i - ^/iV */2 , (9c) 

where W = dx + idy and V* = dx — idy are the com¬ 
plex derivative operators. Some of the notation has been 
changed as compared to the one used in [H to high¬ 
light the comparison with the equations obtained from 
the Boltzmann approach. Specifically, we have tried to 
keep notations as close as possible to that of [ 2 ^, using 
subscript ’S’ or ’B’ for the sets of continuum equations 
derived from the Smoluchowski and Boltzmann equa¬ 
tions respectively. A few further coefficients needed to 
be added, like the coefficient A' in Eq. (Hi), for which we 
kept the original notation since no similar term appear 
in [ 2 ^ (see below). Additionally, we have used Cs = A, 
ng = 3A"/5, k's = A'748, xs = A"/24, and x's = 3A"/5 
where each of the A’s are proportional to the length of 
the rods and the square of the self-propulsion speed. The 
parameter ps = 4iA/{(p/pc — 1) controls the stability of 
the uniform isotropic state, with pc a critical density that 
scales inversely with the square of the length of the rods 
and the square of the self-propulsion speed, Uq, i.e., longer 
rods and faster rods are destabilized at lower densities. 

On the other hand, the equations found in [ 2 ^ using 
the Boltzmann approach are (note that in vq was set 
to ijl 

dtp + voRe(V*fi) = 0, (10a) 

dth = -{a - li\h?)h + CB/r/2 - y (Vp + VV 2 ) 

+ |/2 V/2 , (10b) 

dtf2 = (PB - i\f2?)f2 - y V/i - ^/rV/2 
- ^/2V*/l - ^/lVV2 


The complex fields fi and /2 are related to the vectorial 
polar field P and to the tensorial nematic field Q as 


pP 


'vim/i J 


pQ 


1 

2 


Re /2 Im /2 \ 
Im /2 -Re /2 J ’ 


( 8 ) 


+ ^A /2 -b w/f -b TI/ 1 IV 2 • (10c) 

where A = VV* is the Laplacian. Note that both co¬ 
efficients p 5 and pb are positive at high density and/or 
low noise, leading to a linear instability of the isotropic 


where p = fo denotes the density field. To compare the 
continuum equations obtained in and in |20l. . we 

rewrite them in a common notation. In two dimensions, 
the complex notation /i and /2 is much more convenient 
than the more standard vectorial/tensorial notations. 

The continuum equations obtained in read (ne¬ 
glecting terms arising from spatial diffusion for ease of 


1 In 0 , a ’renormalized’ speed has been found in the equation 
for fi instead of uq, but we neglect this correction at this stage 
to simplify the discussion. We come back to this point in Sec. HY] 
^ As mentionned in [ 2 ^ . Eq. mOcH as published in [20| includes 
a misprint, as the term ^A/2 was erroneously written yV^/2, 
which is incorrect since complex operators are used here instead 
of vectorial ones. 
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state (/i = /2 = 0 ) towards the onset of nematic order 

(/2 ^ 0 ^ 

In spite of a number of similarities, the sets of equa¬ 
tions (I9al9bl9c|) and (llOallObllOcI) do exhibit some differ¬ 
ences. We now highlight them and identify their origin, 
summarizing the detailed analysis given in subsequent 
sections. 

1. Equation (1^ for / 2 , obtained from the Smolu- 
chowski equation, does not contain a saturating 
non-linear term I/ 2 P /2 that is needed to cutoff 
the linear instability and obtain an ordered state 
(although this term was added by hand in m). 
When positional diffusion is neglected, Eq. ([9c]) 
does not contain a diffusion term A /2 either. Both 
these terms are present in Eq. (llOcI) thanks to the 
Ginzburg-Landau closure ansatz used in [ 2 ^ that 
includes higher order modes than the closure used 
in [I 3 II This closure also yields terms I/ 2 P /1 and 
/ 2 V /2 that are included in Eq. (llObp . but not in 
Eq. (l9bl) . 

2. The terms //V/i, /iV*/i, and /iV// are present 
in Eq. (I9bl) , but not in Eq. (llObI) . Similarly, a term 
/ 2 V// is present in Eq. (1^ . but not in Eq. (llOcD . 
These terms arise because the Smoluchowski equa¬ 
tion derived in [13 incorporates excluded volume 
forces and torques arising from the finite size of the 
particles. In contrast, the Boltzmann equation used 
in Ref. considers effective alignment rules be¬ 
tween point-like particles, and thus includes only 
torques (acting during collisions), but no forces. 
We note that the finite size of particles could also be 
incorporated in a Boltzmann approach This 
difference is therefore at the level of the underly¬ 
ing microscopic model and is not associated with 
differences in the closures used. 

3. The terms ff and I/ 1 P /2 are present in Eq. (llOcL 
but not in Eq. (1^ . As we will show below, this 
difference is also due to the different microscopic 
models used by the two sets of authors, which re¬ 
sults in a different symmetry of the collision kernels 
considered in the two kinetic equations. Baskaran 
and Marchetti considered long, thin rods whose col¬ 
lisions are described by a kernel with pure nematic 
symmetry. Peshkov et al. considered instead point 
particles with nematic alignment and an effective 


® In [T^ . the fast modes /a and /4 were discarded. This approach 
is in line with standard practice in kinetic theory of molecular or 
granular gases, where a closure relation is obtained by constrain¬ 
ing the one-particle phase-space distribution to be a function of 
the slowly relaxing fields (the so-called normal form of the distri¬ 
bution) 1 ^ . A higher order closure is, however, needed in both 
passive and active systems to derive the nonlinear terms that 
yield the ordered phase (see | 21 || in the case of active systems). 
Note that this was recognized by Baskaran and Marchetti in llTl 
where the nonlinear term I/2P/2 was added by hand. 


circular excluded volume, resulting in a collision 
kernel that contains both terms of nematic and po¬ 
lar Fourier components. We will see below that 
considering infinitely thin rods in the Boltzmann 
framework also leads to a kernel with nematic sym¬ 
metry. 

In this section, we have summarized the differences 
between the two models both at the kinetic and at the 
hydrodynamic level. We have also briefly identified the 
origin of these differences. The technical aspects of these 
conclusions are unfolded in detail in the subsequent sec¬ 
tions. In section ITTTl we show that a strong formal anal¬ 
ogy emerges between the Smoluchowski and Boltzmann 
equations when considering the limit of point-like parti¬ 
cles with interactions reducing to alignment rules. The 
continuum equations obtained for the order parameters 
have essentially the same structure in both cases, and 
take precisely the same form in the limit of thin rods. 
Then, in section HVl we explain how further terms emerge 
when including forces in addition to alignment torques 
in the Smoluchowski equation, thus obtaining a general¬ 
ization of Eqs. (I9bl9cl) . taking into account the relevant 
(lowest order) nonlinear terms. 


III. POINT-LIKE PARTICLES WITH NEMATIC 
ALIGNMENT INTERACTIONS 

In this section, we derive continuum equations both 
from the Smoluchowski and Boltzmann equations for 
point-like particles interacting only via alignment in¬ 
teractions (i.e., torques). In the case of the Smolu¬ 
chowski equation, this corresponds to only incorporat¬ 
ing the mean torque given by the first term in Eq. ©, 
but neglecting the mean forces given by the second term 
of Eq. ([3]) and Eq. (|T|). In the case of the Boltzmann 
equation only torques were included from the outset in 
Eq. ([2|). By point-like particles, we mean that we con¬ 
sider the limit where explicit excluded volume contribu¬ 
tions are neglected, although the particles still have a 
“shape” that defines the region of interaction and there¬ 
fore determines their collision rate. This shape is chosen 
as needle-like in the work by Baskaran and Marchetti and 
as circular in the work by Peshkov et al.. 


A. Smoluchowski equation 

Following an approach similar to the one used [23 - 
[ 2 ^ for polar self-propelled particles, we first consider the 
force-free Smoluchowski equation, in which interaction 
between particles appears through the average torque r 
exerted by neighboring particles (from now on, we set 

Vo = 1 ) 


dtf + e{e)-Vf = Dnd^ef-de{fT). (II) 
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The local average torque r can generically be expressed 
as 

T(r, e,t)= j dr' J dd'klv' - r, O', 9) fir', 9', t) . (12) 

The gradient term given in Eq. ([3]) is obtained by ex¬ 
panding Eq. (fT^ in gradients and assuming that |r — r'| 
is of order of the length of the rods. To simplify the 
comparison with the Boltzmann approach that consid¬ 
ers point particles, we neglect these terms here and sim¬ 
ply write K{v' — r,d',0) = ^(r' — r)K{9' — 9), where 
we have also taken into account rotational invariance. 
We assume space reversal symmetry, which leads to 
K{—9) = —K{9). We also assume that K{9) obeys 
a nematic symmetry K{9 tt) = K{9), which is valid 
for thin rods undergoing collisions, when cap-on-cap col¬ 
lisions can be neglected El- Hence the local average 
torque reads 

r(r, 9,t)= J d9'K{9' - 9)f{r, 9', t) (13) 

Passing to Eourier components, we obtain the following 
hierarchy of equations: 

dtfk + + iv7fe+i (14) 

., oo 

= -^kfk + ^ E 

q— — co 

where = Duk, and 

kg= f e^<i^K{9)d9 (15) 

J —TT 

Due to the nematic symmetry of the interaction, all odd 
Fourier modes K 2 m+i = 0. Defining 

Jig = ^k.g (16) 


distribution, that is, for small values of the Fourier har¬ 
monics fk {k > 1), Eq. (1171) can be rewritten to linear 
order as 

dtfk + Ivfk-i + iv7fe+i = 9k fk (19) 

with fil = (J^j, +Jlo)P~kJk- density, the linear 

coefficient is negative. Whether /xf becomes positive 
or not at higher density depends on the sign of Jlk + ^ko- 
The precise values of these coefficients depend on the de¬ 
tails of the chosen interactions, which we do not specify 
explicitly here. We however assume that interactions fa¬ 
vor nematic alignment, resulting in 7^2 + ^ 2.0 > 9. As 
a result, /xf becomes positive above a transition density 
Pt -or equivalently, below a given threshold value of Df 
at fixed density. For p just slightly above pt , /xf is posi¬ 
tive and small, leading to a slow instability of the state 
/2 = 0, a regime in which the dynamics of the system can 
be reduced to that of a few coupled modes. In constrast, 
we assume that J^i + j£o < 0 since interactions do not 
favor polar order, so that /xf remains negative. 

Following [^[^, we introduce a truncation procedure 
close to the instability threshold of the linearized equa¬ 
tion. We use the following scaling ansatz, with e a small 
parameter: 

P-Po ~/i ~/2 ~ e, f 2 k-i^f 2 k^e^ (fc > 0) (20) 

To lighten notations, we further introduce the coefficient 
Clg defined as 

Cfe.g = Jk,q + Jk,k-q ( 21 ) 

The nematic symmetry, which implies Jig = 0 for odd 
g, in turn leads to Cf ^ = 0 if A: is even and q is odd: 
^2m,2l+l ~ 9- 

After truncation of Eq. cm to order for fc = 1 and 
fc = 2 , one obtains the equations governing the evolution 
of /i and / 2 , 


we get for the Fourier transform of the Smoluchowski 
equation 

dtfk + iv/fe_i + iv7fc+i (17) 

OO 

= ~^kfk + E/ dlgfk-qfq 
q— — cio 

Note that due to the symmetry K{—9) = —K{9), Jig is 
real. In addition, jf ^ = 0 for odd q due to the nematic 
symmetry of the interaction. 

For fc = 0, Eq. (ITTl) directly leads to the continuity 
equation 


atp + Re(V 7 i) = 0 (18) 

To derive closed equations for fi and 7, we need to re¬ 
sort to an approximation scheme. For an almost isotropic 


9t/i + ^(V72 + Vp) =/xf/i (22a) 

+ Cf 2 /i 72 + Cf 3/273 

dtf2 + i(V73 + V/i) = plh + ClJih (22b) 


These equations are not closed, as they also involve the 
higher order harmonics /a and 7 . We thus make use of 
Eq. (1171) for fc = 3 and fc = 4, truncating them to order 
since /a and 7 appear in Eqs. (I22a|) and (I22bl) only 
in space derivatives or multiplied by another small field. 
We then obtain 


h 

h 


1 <^ 3^2 

:^V /2 - -^hf2 

2/xf /X* 


CI 2 

' 2 pf 


5/2 


(23a) 

(23b) 








6 


Injecting Eqs. (j23al) and (I23bl) in Eqs. (I22al) and (j22bl) . 
we obtain closed equations for fi and / 2 , 


dtfi = Mi/i - 


r'S nS 

S £ ^1,3 '-^3,2 I £ |2 


1 / 2 IV 1 + cyu2 


Ms 






/ 2 V /2 


at/2 = Mf/2 - ^^l/2p/2 - iv/i 


2mI 


- A^/2 

4M3 


cio 


(24a) 


(24b) 


B. Boltzmann equation 

In the case of the Boltzmann equation, interactions are 
described as instantaneous collisions, and the equation 
reads 


at/ + e(0). V/ = /dif[/] + /coi[/] (25) 

where Ida and Icoi are defined in Eqs. m and ([6|) respec¬ 
tively. The dynamics is defined such that the angles 0[ 
and 02 after collisions are given by 

e[ =-^{91,62)+ m, 02 = 4 '( 02 , 0 i) + r72 (26) 

where 771 and 772 are independent noises drawn from P{ri). 
We assume here that 'l'(0i, ^ 2 ) favors nematic alignment. 
Using symmetry prop erties. 47(01,02) can generically be 
parameterized as 

4'(0i, 02) = 0i+i?(02-0i) (27) 


where H is an arbitrary function which encodes the 
nematic symmetry of the interaction, and is thus tt- 
periodic. After expansion in angular Fourier series, one 
finds for the Boltzmann equation Sim, 

dtfk + + iv7fe+i (28) 

00 

= -D^h+ E 

q— — cio 

where the coefficients and are given by 


with 


Uf = -Afl(l-Pfe) (29a) 

•^k,q = Pklk,q — lk,0 (29b) 


Pk = 

^k,q — 


d7]P{r]) 

o 

deKeiO) 


(30a) 

(30b) 


Note that Eq. (l28l) is formally identical to Eq. (flTl) . The 
only difference, apart from the specific values of the co¬ 
efficients, is that is a priori nonzero for all {k,q), 
as it is not constrained by a nematic symmetry. In the 
explicit model considered in where particles are con¬ 
sidered to have an almost circular shape of diameter do 
(the interaction radius), one has (setting uq = 1) 


Kb{9' 


0) = 4do 



(31) 


which is not invariant under the nematic symmetry 0 
9 + TT. This absence of nematic symmetry is due to the 
fact that K{9' — 9) corresponds to the collision rate of 
polar particles, which depends on the velocity difference 
of the particles. However, for non-circular particles, a 
simple heuristic generalization of Eq. m can be pro¬ 
posed, modulating the collision rate by an orientation- 
dependent interaction radius d{9' — 9). In the limit of 
infinitely thin rods, one finds (in the frame moving with 
the particle of orientation 0), 


d{9' 


9)=do 



(32) 


so that the collision kernel obtained in this limit. 


K 


thin 

B 


{9' 


9) = 4d(0' - 0) 



2do| sin(0' — 0)| 


(33) 

indeed obeys a nematic symmetry. One thus recovers, in 
the limit of infinitely thin rods, the property 21+1 ~ 
0, as in the Smoluchowski case. 

The derivation of the continuum equations for the po¬ 
lar and nematic order parameters follows exactly the 
same lines as in Sect. IIII A[ and one finds (see for 
details) 


dtfi = /if/i - 


/^B 

^1,3 ^3,2 


B l/2|7l+Cf2/l72 
/if 


-ivp-iv72 


/^B 


f^B ^B 

o i? B £ 2,4 4,2 \ £ \‘2 £ 

dtf2 — f^2 h --1/2 I /2 


^B /^B -I 

/is ^ 


I 4 ^B j*2 '-^2,S'-^S,2 

+ 2^2,iJl ~ 


1 


(34a) 


(34b) 


4/i^ 


■A/2 


^'1v*(/i/2) + ^/rv/2 


B 


2/if 


2/if 


These equations have essentially the same form as 
Eqs. (I24all and (I24bll . The additional terms ff, I/ 1 P /2 
and / 1 V /2 that appear in the equation for 7 are due to 
the lack of nematic symmetry of the kernel Kb{9' — 0). 
These terms vanish in the limit of infinitely thin rods, 
where one recovers the same equations as obtained by 
Baskaran and Marchetti on the basis of the Smoluchowski 
equation. In other words, when only torques, but no 
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forces are included in the kinetic equation, the two ki¬ 
netic equations yield the same continuum equations and 
the differences in the published continuum equations ob¬ 
tained by the two approaches are entirely due to differ¬ 
ences in the microscopic models. 


IV. PARTICLES WITH INTERACTION 
FORCES AND TORQUES 

To understand the remaining differences between the 
two set of published continuum equations we need to ex¬ 
amine the remaining terms in the Smoluchowski equa¬ 
tion. Specifically, Eq. (l9bl) contains terms coupling the 
polar order parameter with its space derivative that are 
not in Eq. (llObI) . In terms of the complex notations used 
here, such terms yield a linear combination of /* V/i, 
/iV*/i and fiVf*. Also, the term / 2 V /1 appearing in 
Eq. (1^ was not obtained in Eq. (I24bl) . 

These additional terms arise because the Smoluchowski 
equation derived by Baskaran and Marchetti contains two 
additional terms as compared to the Boltzmann equation 
used by Peshkov et ai. The first one is the mean force 
given in Eq. o that describes momentum transfer in 
a collision. The second one is the second term in the 
torque given in Eq. ©, that arises from the difference 
in the position of the center of two colliding rods, and 
hence incorporates the finite size of the particles. Note 
that both terms could also be incorporated in a Boltz¬ 
mann approach. Both terms are of first order in the spa¬ 
tial gradients and yield terms of the same symmetry in 
the continuum equations. In this section we discuss the 
terms arising form the mean force, while those due to 
the nonlocal torque are discussed in Appendix In ad¬ 
dition, Baskaran and Marchetti incorporated positional 
diffusion in [13, leading to the following Smoluchowski 
equation, with the short notation e = e{9), 

dtf + Cadaf = DoAf + Di 

+ Dndy-de{fT)-V-{fF) (35) 

where Dq and Di are the isotropic and anisotropic dif¬ 
fusion coefficients, and F and r are the force and torque 
given in Eqs. © and ©, respectively (but we will ignore 
here the second term on the right hand side of Eq. 
and discuss it in Appendix lAT) . 

Starting from the Smoluchowski equation (1551) . the fol¬ 
lowing continuum equations (rewritten here in the cur¬ 


rent complex notation) have been obtained [13, 
dtp + Re(V*/i) = DpAp + :^Re(V*V 2 ) (36a) 

dtfi = -Dnh - \^P - ^VV2 + A/r/2 
-A'(/rV/i + /iVVi-/iV/r) 

+ i(A + Dsj,i)Ah + \{Db- Dspi)V‘^fl (36b) 

dtf2 = Psf2 - y V/i + ^V^p 

o\// \'f 

- ^(/rv/2 + /iv*/2) - ^(2/2 VVi+/ 2v/r) 

(36c) 

where vq is a renormalized speed. In the following, we 
aim at rederiving the generic form of Eqs. (l36al36bl36cL 
including the relevant nonlinear terms. 

On general grounds, the kernel G{9,9') in Eq. © as¬ 
sociated to the force exerted by a particle of orientation 
9' on a particle of orientation 9 can be decomposed onto 
the directions parallel and perpendicular to e{9), 

G{9,9') = G||(0' - 9) e{9) + G^{9' - 9) e^{9) (37) 

with ej_(0) = e{9 + ^). The scalar functions G|| and 
G_l depend only on the angle difference, by rotational 
symmetry. For non-chiral particles, the force G(0, 9') 
obeys a reflection symmetry, characterized by 


G||(-0) = G||(0), G^{-9) = -Gde) (38) 

The angular Fourier transform of fF can be decomposed 
into parallel and transverse contributions: 


J 

d9 9, t)F{r, 0, t) = f'' + F,^ 

— 7T 

(39) 

with 


pTT 


f]} = 

1 d9e^^^f{9) j 

l-K 

1 d9'f{9')G\\{9' -9)e{9) 

— TT 

pTT 

(40a) 

II 

J d9e^^^f{9) J 

i (i07(0')Gj_(0'-0)e_L(0)(4Ob) 

— TT 


where we have written simply f{9) instead of /(r, 9, t) to 
lighten notations. After some algebra, we can eventually 
write V • [fF) in the form 

00 

V-{fF)= Y. 

q— — oo 

CO 

+ ^ M,V(/,A-_1_,) (41) 

q— — co 

where 

^9 = ^(Gi'+*G^) (42) 








with Gq and the Fourier transforms of G|| and G± 


■TT 



1 d 6 »e*'=®G||( 6 ») 

— TT 
pTT 

(43a) 


1 d 6 )e"G_L( 6 ») 

— TT 

(43b) 

1 


The Fourier coefficients Gq and G^ are constrained by 

symmetries. The coefficient Gq is real and is purely 
imaginary due to the reflection symmetry Eq. (1381) . It 
thus follows that Mq is real. Note also that = Gq 

and Gfq = —G^, a property that we have used in 
Eq. (HTI) . 

In Fourier transform, the Smoluchowski equation thus 
reads (see m for a similar derivation of the anisotropic 
diffusion terms) 


1 1 D ^ 

at/fe + -V/fc_l + -VVfc+l=^oA/fe + ^(VVfe- 2 +V*Vfe+ 2 )-^f/fe+ ^ Jlqfk-qfq 

q——oo 

+ M_qV*{fqfk+l-q) + ^ MqV{fqh-l-q) (44) 


As we show below, applying the closure scheme used in 
Ref. [13 one obtains continuum equations for p, fi and /2 
that have the same form as those of [13. Of course, the 
detailed expression of the coefficients differs since we have 
neglected the gradient term in the torque, but we were 


interested here only in the structure of the continuum 
equations. 

To examine the effect of the closure, we truncate 
Eq. (l44l) to order for p, fi and / 2 , yielding 


atp + Re(V*/i) =44oAp+^Re(V*2/2) + 2(Mo + M_i)Re(V*(p/i)) +2(Mi + M_2)Re(V*(/i72)) (45) 

dtfi + \v*f 2 + (i - Mop)Vp = pf A + GoA/i + ^VVr + Cy;f 2 + Gf,3/273 

+ (Mo + M_ 2 )VAp/ 2 ) + (M 2 + M_2)V|/2|2 + M_iV7i" + (Ml + M_i)V|/ip (46) 
dtf2 + i(V73 + V/i) = plf2 + G 0 A /2 + ^V^P + GlJ^U + (M_l + M_2)V*(/72) 

+ (M _1 + M2)V(/i72) + (Mo + M_3)VAp/3) + (Mo + Mi)V(p/i) (47) 


For /s, Eq. (l23aD is changed into 

/3 = 4 f T - (^0 + ^2)p) V/2 - %/i/2 (48) 

M 3 \ ^ / M 3 

while one recovers Eq. (I23bll for A. Using these last 
equations as closure relations, we find for fi and for A 

J 


dtfi = Mf/i - - Mop) Vp - iv72 + GoA/i + ’' iM'/i + CIJU 2 + 7/2 V /2 

+ 2 A/_i/iV 7 i + (Ml + M_i)(/iV/i* + /rv/i) + (Mo + A/_2)V*(p/2) + (M2 + M_2)f2^f; ( 49 a) 

dtf2 = 4f2 - ^^0^\f2?f2 - iv/i + (Mo + M,)V{ph) + GA /2 + 

+ x(/iV 72 + / 2 V 71 ) + (M_l + M 2 )(/i*V /2 + / 2 V/ 1 *) 


(49b) 
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where the coefhcients are given by 

7 = ~ ("^0 + M 2 )p^ + {M 2 + M- 2 ) (50a) 

D = Dq - (2 ~ -^- 3 )^) (2 ~ 

(50b) 

X = (2 ~ -^- 3 )^) + {M-i + M_2) (50c) 

Note that, unlike the torque, the force kernel used by 
Baskaran and Marchetti [l^ does not have the nematic 
symmetry, as in general G{6+tt, O') ^ G(0, O'). Were this 
symmetry present the Fourier coefficients would satisfy 
^ 2 n = ^ 2 n = leading to M 2 „ = 0 for any integer n. 


V. DISCUSSION 

In this paper, we have examined the differences be¬ 
tween the continuum equations for interacting self- 
propelled rods previously obtained in the literature from 
the Smoluchowski equation and from the Boltzmann 
equation in The differences, not surprisingly, arise 
either from (i) the use of different microscopic models 
or (ii) the use of different closures of the moment ex¬ 
pansion of the kinetic equation. Concerning the model: 
Baskaran and Marchetti considered long, thin rods with 
finite and anisotropic excluded volume (hence incorpo¬ 
rating momentum transfer in a collision and nonlocality 
on the scale of the difference in position of the colliding 
rods), while Peshkov et al. considered point-like particles 
with nematic alignment rules and circular interaction ar¬ 
eas. Additionally, while Baskaran and Marchetti used a 
simple truncation that neglects all moments higher than 
the second, Peshkov et al. employed a more sophisticated 
closure that allows one to derive in particular the nonlin¬ 
ear term responsible for the onset of the ordered state. 
We show here that when the same microscopic model and 
the same closure are used, the Smoluchowski and Boltz¬ 
mann approach yield the same continuum equations, al¬ 
beit with different microscopic expression for the param¬ 
eters. This is perhaps not surprising, but it is reassuring 
to demonstrate the equivalence for these nonequilibrium 
systems. 

One interesting result is that even in the simplest case 
of point-like particles with nematic aligning interactions, 
where only strictly local aligning torques are included, 
the shape of the interaction region controls the symme¬ 
try of the collision frequency and the symmetry of the 
interaction kernel. The needle-like particles considered 
by Baskaran and Marchetti yield a collision kernel that 
contains only even Fourier components in the angle de¬ 
scribing the difference in orientation of the interacting 
particles, and hence has pure nematic symmetry. The 
circular point-particles considered by Peshkov et al. yield 
a collision kernel that contains both even and odd Fourier 


components, and hence has mixed polar and nematic 
symmetry. This leads to additional terms in the equa¬ 
tions of Ref. [13 not obtained in Ref. [l3,13 • We ex¬ 
pect that these additional terms will generally be present 
when one considers finite-thickness rods, so that cap-on- 
cap collisions are not negligible. Numerical integration 
however indicates that at low density, these terms do not 
play an important role in the behavior of the equation. 

The terms /fV/i, /iV*/i and /iV/j" present in the 
work by Baskaran and Marchetti, but not in that of 
Peshkov et al., arise from both the mean force given in 
Eq. ([4]) and the nonlocal contribution to the torque in 
Eq. ([3|. Here the discussion of the latter has been rele¬ 
gated to Appendix A, but both contributions to the ki¬ 
netic equation must be included to consistently evaluate 
the coefficients of these terms, as shown in [i7|. These 
terms encode the fact that the polarization in the hard 
rod models is actually the physical flow velocity of the 
fluid of rods. They include the convection term P • VP 
that is the analog of the convective nonlinearity in the 
Navier-Stokes equation of passive fluids and a V|Pp con¬ 
tribution to the pressure that arises from self propulsion. 
The mean force and nonlocal torque also yield additional 
terms Re(V*(p/i)) and Re(V*(/f/ 2 )) in the continuity 
equation (H 51 ) . These do not appear in Eq. (ISal) because 
nonlinear terms in the hydrodynamic variables in the 
continuity equation were neglected in [la. Isotropic and 
anisotropic diffusion terms, as appearing in Eqs. (I36a|) . 
(I36bl) and (l36cD , are also reproduced when positional dif¬ 
fusion is taken into account, see Eqs. (1^ . (|49a|) and 
(|49bp . Finally, a renormalization of the velocity vq = 1, 
in agreement with Eq. (I36cll . is obtained in Eq. (I49bp . by 
expanding the term V(p/i). 

Note that the non-locality of the interactions, only con¬ 
tributes to order to the average force. These terms 
are therefore neglected in the truncation procedure. The 
only terms of order in the equation for /fc (fc = 0, 1, 2) 
are of the form pAfk, pV^ fk -2 and pV*^ fk +2 (for k = 0 
only in this last case) and thus would just renormalize 
the existing diffusion terms. 

Finally, let us mention that the non-locality of inter¬ 
actions (or, in other words, the finite excluded volume 
of particles) can also be accounted for in the Boltzmann 
framework. The corresponding study goes beyond the 
scope of the present note, and will be the subject of a 
future publication [^ . 
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Appendix A: Non-local average torque 

In this appendix, we briefly sketch the derivation of 
the Fourier transform of the Smoluchowski equation in 
the case of weakly non-local interactions affecting the 
torque only. We start from the non-local average torque 
Eq. (fT^ , and first rewrite the kernel using the change of 
variables r' — r = se (</>). From rotational symmetry, we 
obtain 

K{r' -r,e',e)=K{s,9' -e,(j)-0) (Al) 

Expanding the torque to first order in gradient, we get 

T{e)= f de'K,(e'-0)f{9') 

J —TT 

/ TT nTT 

d9' / d(j)K2{0' -9,(1)-9) e((^) • V/( 6 »') 

-TT J —TT 

= ti{9)+T2{9) (A2) 

where we have defined 


/ TT pOO 

d(l) / ds sK{sj6'— 9^(j) — 9) (A3a) 

-TT Jo 
pCO 

K2{9'-9,(j)-9) = dss^K{s,9'-9,(j)-9) (A3b) 

Jo 

Note that Eq. is recovered by introducing 

K2{9',9)= f d(l)K2{9'-9,(k-9)e{(^) . (A4) 

J —TT 


From Sect. IIII Al we know that the Fourier transform of 
Ti (9) is given by 

/ TT 

d9e^^\,{9)= E (A5) 

■TT ^ 


A similar calculation for the Fourier transform of T 2 yields 
f d9e^'^^T2{9) = E (A-,/fc+i-,V*/, - LJk-i-gVfg) 

J-TT q 

with 


(A 6 ) 


Ln = K„ I 

^ 47r ^ 


(A7) 


and 


f d9i f (A 8 ) 

J —TT J —TT 

Note that due to the nematic symmetry K2{9i +t(,92 + 
tt) = K2{9i, 6 * 2 ), one has L 2 n = 0 , hence the sum may be 
carried over odd q’s only in Eq. (IA 6 I) . Pluging the ex¬ 
pression (IA 6 II of the torque into the Fourier transform of 
the Smoluchowski equation (with force term) eventually 
leads to 

c^tfk + + 2'^*fk+i 

= DoAfk + ^iV^fk-2 + V*2/fc+2) 

00 

-Dlfk+ E 

q— — cio 

+ + M-q + Mq-k-l)fk+l-q^* fq 

Q 

+ Y,{kLq + Mq + Mk-l-q)fk-l-qVfq (A9) 
<? 

The resulting equation thus has the same structure as 
Eq. (|44)l . after expansion of the derivatives of products 
and relabelling. Only the values of the coefficients differ. 
This however leads to a reweighting of the different terms 
which generically breaks the symmetry between terms 
like / 1 VV 2 and / 2 V*/!, or /*V /2 and / 2 V/* observed 
in Eq. (I49bp . 


q— — oo 
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